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1. Introduction 

Exchange rate fluctuation, values of shares, fiitures and options are random or quasirandom processes, 
and many "objective" (or perceiving as such) and random or seeming (perceiving) as random factors influence 
on them. 

Here a number of questions arise. Whether it is possible on the basis of analysis of limited numerical 
sequence of financial instrument cost to draw a conclusion about the nature of statistical processes underlying 
the observed event, specifically, to decide a question about stationarity (quasi-stationarity) of the process, to 
restore probability density upon observed data, to calculate probability of certain events etc.? 

Preceding the result of suggested work, it is necessary to point out that the answer is positive, specifically, 
within the framework of not so critical axiomatics, creation of desired statistics on the basis of analysis of 
random or quasirandom sequences of financial instrument costs is possible. 

Furthermore, according to the author, solution of the problem is essentially eased because we do not 
speak about simulation of complicated natural processes, but only about simulation of decision making method 
used by market participants, who have the same initial information and comparable computational and 
analytical capabilities, which pursue well-defined and, generally speaking, obvious target. 

2. Model of nonlinear stochastic wavelets (NSW model) 

Let random process X{t, CO) be the cost of shares, options and futures in the point of time t (^ R , 
where CtX^Cl is a point in probability space with measure function P . Let's introduce a population of 
analyzing functions-wavelets y^i^{t) = 2 ' y/{2^'t - r) , meeting the condition i//{t) t:^l} r\l} , y/{t)dt = 

and forming the orthonormal basis for L {R) with compact carrier. Small value of index j conforms to low- 
frequency component of the random function Jf(?,<Z>) , and the large one to high-frequency component [1]. 

Let's express the random function in the form of linear combination of functions 
Yi(t)=\X(t,af)if/i^(t)dt. We designate Y(t) = (Y^(t),..YJ(ty) , where - / is maximum order of wavelet 



transformation. We presume that new J - dimensional random vector Y(?) is in accord with stochastic 
equation of Ito: 

dY{T) = ¥{Y{T))dT + G{Y{T))-d&, (1) 

where dcOi{t) is random Wiener process <da)^{t)>=0, <da)-{t)da>j{ty)>= 5y5{t -t^), i,j = l...J. 

Symbol < . . . > means averaging as per random process X(t, of) with measure function P and 
F(Y(r)), G(Y(r)) are drift and diffused components accordingly, we assume thatG,^(Y(r)) = SyG(Y{T)) . 

Conventional transition probability /(y(0/y(^o)) ^^^ be represented as a path integral [2]. Nonlinear 
effects can be taken into account using diagrammatic method. Specifically, we can say that for nonlinear model 
the spectral potency (Fourier transform of two-point twin correlation fiinction), some times has power 
asymptotics, i.e. observed "fractalness" of signal is only the result of nonlinearity of the model. 

In single-mode case J = 1, time-independent density of distribution defined as 
P(Yi < yi) - \ fs(yi)^yi ^^^ t)^ represented in the form of quadratures: /^(ji)*^ , where 
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In general J - dimensional case fulfillment of additional requirement, i.e. condition of detailed balance 
[3], is required and it is natural. However, if stationary solution exists, then it always can be represented in 
exponential form/, (y) we ^'^'^''^^ where ^i(y)-"^4(y) are certain functions of random variables y. This 
result can be received based on the principle of maximum informational entropy [3] or using inequation of Rao- 
Cramer. In this case the moments < ^/'■i(y)"-V'^A:(y) ^ ^^^ limitations (links), under which we find extremum 

<ln/.(y)>. 

Let Hf,{x^j) be adequate set of orthogonal functions, for example, Ermiht polynoms of k order. 
Orthogonal decomposition F(Y) , G(Y) per polynoms Hi^(Y^) is given by [2]: 

F(Y) = XZ^::f^.(^«)-^/(^p) 

n...p k...l 

GA^) = YLQ::ijH,{yn)-H,{y,) 

n...p k...l 

(2) 



Aim of further stochastic simulation is reconstruction of motion equation (1) and conditional probabilities 
within the framework of made suppositions based upon available observations (experimental information) or, 
that is equivalent, determination of decomposition values in proportion (2). 

To solve the problem we shall use extremal principle, i.e. we shall define the value /t-^ j '' , Ql'"":'' in 

such a manner that minimize functional that are norms of differences of observed and simulated values 
averaged per random process with measure function P . After differentiating these expressions with respect to 

^Lj''' Ql "Y ^ multiplying by Hi^(y^) and integrating on y„...yp we have two systems of linear equations 
for unknown variables yil^j'' , Qk'""/ ■ 

Solution of these systems gives us expression for drift and diffused factors of non-Markov stochastic 
equation, i.e. the most probable path function of the process and statistic of fluctuations round it. We shall name 
corresponding sequence of random variables as setting up sequence of T^ length. 

On the basis of these results we can construct optimal policy of portfolio management. It follows from 
compactedness and condition y/{t)dt = that factors of wavelet decomposition characterize mean change of 

value X(t, co) on the length of carrier 2 "^ Tq . 

Therefore, in elementary case, moment of optimal purchase of securities can be represented as 
t = vaf{-dy^>0,P^>\-a^) and moment of optimal sale as ^ = rnf(-(i_Vj < 0,/^. < «[), where 



P^= /, {y)dy . Sign of dy^ has been chosen for Haar wavelet, and value is risk level. The model under 

—00 

consideration has non-Markov nature. 

Conditions P^ >\-a^,P^ <a^ are required for exception (on the average) of the purchase in case of 
7i (r„_i ) < 7i (r„ ) , - J7„ < and sale in case of Y, {t„_, )>Y,{t„),- dY„ > . 

-HX> 

Alternative variant is use of convolution f^ (z) - f^ (y^ )f^ j. (y^ + z)dy^ . In this case criteria 
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P^ > 1 — «, , P, < «j mean resale- or repurchase moments. 

In both cases it is necessary to exclude statistically nonsignificant differences between stationary 
distributions in time displacement t ^^t + T . 

The last is possible when using non-parametric Kolmogorov criterion with risk level a 2, i.e. synthesized 
stationary distributions fXy\ )■> fXy\ ) in displacement t ^^■t + T with probability /? = 1 - QTj are considered 



to be indistinguishable, if tnax\K(y^^(t^) — K'(yi(t„)^< k(a2)/ vN , where K,K' are cumulative 
probabilities for this condition and condition received as a result of displacement T . 



The constant k(a2 ) ~ 1 depends on risk level ttj , A^ is the number of experimental points (samples) per 
which approximation of stationary distribution has been constructed. Conditions of purchase/ sale of financial 
instruments in NSW model are shown schematically in figure 1. 



X(t) 




Ti 



Figure 1. Conditions of purchase/ sale of financial instruments in NSW model; here z^is compactedness 

interval of analyzing wavelet, T^ is length of setting up sequence, Tis interval of quasi-stationarity of 

distribution with risk level «j , Tj is interval of parcel redistribution decision making, and A, B are points of 
purchase and sale of financial instruments. 

3. Experimental proof of NSW model 

Data of three types that are stock prices of Russian emitters (Sberbank, SurgutNefteGas and Gasprom) 
and stocks of U.S. companies (AT&T Inc, Dell Inc. and the Bank of America) over a period since January 1, 
2009 till December 3 1 , 2009 were used as initial database for model run. 



It was analyzed the non-Markov dual-mode model ( J = 2) , scanning was implemented in increments of 
1 min., calibration sequence amounts to 32-64 indications, and the order of Hermite polynom not to exceed 
K = 3. 

Analyzing function - Haar wavelets ^^ , Daubechies wavelets ^2 » ^3 ^^^ Bettle-Lamerie spline wavelets 
5j , ^2 5 ^3 [1]- Although the last ones do not have compact carrier, they come down quickly enough under large 
time. 

Cost of securities at the start time was assumed asC(^Q) = 1. Objective function of the model is profit at 
the point of time?: Z{t) = C{t)/ €{1^) , decision making for purchase/ sale in accordance with proposed 
criterion and calculation algorithm for values F(Y) , G(Y) at the point of time t > t^ was implemented in 
automatic mode based upon prior information. Risk levels were taken as a^ = aj = 0.05 . 

Results of calculation for three groups of active securities are given in figure 1 . The results are designated 
at intervals of five acts of purchase and sale. 





Figure 2 Profitability when using NSW algorithm within the period since January 1 , 2009 till December 
31, 2009 for stocks of Sberbank (SBER), Gazprom (GP) and SurgutNefteGaz (SNG) (a) and for AT&T, Dell 
Inc. and the Bank of America (b); and risk levels are a^ = aj = 0.05 . 



Average fraction of decisions made (number of purchase and sale acts towards length of purchase and 
sale field) were in the interval 0.01. ..0.025. Daubechies wavelets were optimal for the stocks of Russian emitters 
and stocks of USA companies. 



3. Optimization of security portfolio on the basis of NSW algorithm 

Let z . (t) be the cost of instrument i = l..M at the point of time t in accordance with abovementioned 
decision algorithm. Let x, (t) = \nz-(t + Tq)/ z. (t) be logarithmic earning power and let's define average value 
X. and two-point correlation function at the point of time t > T^: 



1 ' 
^,-- \xXt)dt 

^l l-T, 



A,j^yi(x,(t)-X,)(xj(t)-X^)dt 



Lower estimate of interval Tj is maximum value per / of distribution quasi-stationarity with 
predetermined risk level «[ (fig. 1). Nondiagonal elements A-j take account of intercorrelation of profitability 
of instruments /, j in accordance with NSW algorithm. 



If «;(0 is fraction of instrument / in parcel at point of time t , then in Gaussian approximation average 

M 

profitability of the parcel is Z(t) = 2_jn^(t)X.(t) and its expected mean square is 



M M 



y\t) = Y,T.''.(t)n,(t)x,(t)x^(t). 



Moments of higher order A^ can be set aside, as far as by estimate for all instruments under 



consideration the following is right: 



1, 



I... J 
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Al..A]j 



0.002...0.004«1. 



Let's define the most effective parcel at point of time ^, as a set of instruments i with fraction n^(t) , 
which maximizes value of probability 



P(0) = J i//^(z, y, s)ds -^ extr under limitations ^ «, < 1 . 



eZ{e) 



The parameter 6 >0 implements the compromise between maximum earning power and minimum risk. 
When ^ — > maximum of earning power is preferable, and when 6 ^>\ risk minimum is preferable. Thus, 
the task about the most effective portfolio is the typical task of nonlinear optimization. Necessary and sufficient 



condition of extremum existence is provided with Kun-Takker theorem. To solve the problem we shall use 
gradient procedure. It is possible to take n. =\l M as initial points. 

In Figure 3 there is relation between the average profit and the share fraction in parcel against the time for 
optimal portfolio of AT&T Inc, Dell Inc. and the Bank of America over a period January 1, 2009 to June 30, 
2009. We use Daubechies wavelets. As opposed to Figure 2, cost of the parcel is represented at point of time, 
when decision about redistribution of instruments in the parcel is made. 

Trading was carried out in computer-aided mode without readjustment of the algorithm over a period of a 
half year. Only risk levels a^ , CiTj and compromise parameter 6 between profitability and minimum risk were 
pre -set as the input. 
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Figure 3. a) relation between average profitability of parcel of AT&T Inc, Dell Inc. and the Bank of 
America and time. Risk levels a^ = Ofj =0.05, and compromise parameter 6 = \IA; b) share fraction in 
optimal parcel depending on time. 

In ordinary moment the parcel is shared equally between AT&T Inc, Dell Inc and the Bank of America. 
However, in the process of trading according to NSW algorithm, redistribution of the parcel in favor of the Bank 
of America for the account of Dell Inc. takes place. Accumulated profit of the parcel is less than in case of 
dealing in the Bank of America shares (Figure 2), but it is achieved with the less risk. 



4. Discussion about the results 



In table 1 the NSW algorithm is compared with the following well-known algorithms: Price Channel, 
Bollinger Bands (BB), Moving Average Convergence/Divergence (MACD) and Relative Strength Index over a 
period of January 1, 2009 to December 31, 2009. 



Parameters of standard algorithm were chosen to the effect that maximum earning power is provided. In 
this context, when comparing the algorithms, the well-known models had the advantage a priori; i.e. the 
parameters were chosen according to the process, which had been already realized (in reality it is natural that 
"future" is unknown). 

Proposed NSW algorithm was implemented as close as possible to reality; predictions and control were 
fulfilled on the basic of knowing the previous system status only. Advantages of proposed algorithm are 
essentially for all financial instruments. 

Table 1 Comparison of profitability of various algorithms, January 1, 2009 to December 31, 2009 





PC 


BB 


MACD 


RSI 


NSW 


The Bank of 
America 


1.596 


1.631 


1.273 


1.683 


1.852 


Dell Inc 


1.131 


1.377 


1.311 


1.185 


1.524 


AT&T Inc 


1.433 


1.431 


1.254 


1.439 


1.721 



In table 2, comparison of algorithms was carried out in poor for the Bank of America market conditions. 
There was significant negative trend with strong volatility. Standard algorithms turn out to be unprofitable. 
NSW algorithm allows receiving of small, but stable profitability. 

Table 2. Comparison of profitability of various algorithms for significant negative trends with great 
volatility. The Bank of America, May 1, 2010 to June 30, 2010. 





PC 


BB 


cc 


OC 


NSW 


The Bank of 
America 


0.99 


0.908 


0.918 


0.872 


1.078 



The results analogous to those represented in tables 1 and 2 were received in automatic trading with other 
financial instruments, particularly with futures of Brent (BR-6.10, BR-7.10, BR-8.10, BR-9.10). 



Thus, non-Markov generalization of Ito equation in wavelet-space enables implementation of "short" 
predictions of the most probable path of price change and statistics in the neighborhood of this path. Jointly with 
stationary distribution it enables to optimize the system of automatic trading. 
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